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ABSTRACT 


Recent improvements in rare earth nagnets have made it 
possible to construct strong, lightweight, high horsepower 
DC motors. This has cccasioned a reassessment of electronme- 
chanical actuators as alternatives to comparable pneuratic 
and hydraulic systems for use in flight control actuators 
for tactical missiles. This thesis cevelops a lecw-crder 
Mathematical model fer the Simulation and analysis of brush- 
less DC motor performance. The model is implemented in CSP 
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Electronic commutation based on Hall effect senscr _ frosi- 
tional feedback is simulated. Steady state motor Lehavior 
is studied under both constant and variable air gap flux 
cond tions. The variable flux takes two different foraus. 
In the first case, the flux is varied as a simple sinusoid. 
In the second case, the flux is varied as the sum of a sinu- 


soid and one of its harmonics. 
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I. INTRODUCTION 

Direct current electric motors are jaining wider use in 
space aprilications where long life, hign reliability, high 
torgue and light weight are critical factors. In recent. 
years, advances in magnet technology and materials, esfe- 
Clally the use of rare earth magnets, have made fossible 
Significant imprcevements in the torgue td inertia ratios of 
permanent magnet (PM) motors. Samarium cobalt is cne such 
rare €arth magnetic material. The great potential of the 
rare cé¢arth magnets derive from their inherent high flux 
density and high coercivity properties. ag hers f dl ux mdsens 1— 
ties mean greater developed motor torgue while high coer- 
Civity means greater innate resistance to demagnetization 
which permits thinner magnet sections. tietewOorLactors 1. 
combination result in increased mechanical power and energy 
product with decreased physical size and weight. 

Brushless dc motors belong to this class of permanent 
Magnet motors and thus enjoy their improved torque to size 
enaracteristics. In fact, brushless de motors are signifi- 
cantly smaller than either their ac counterparts or conven- 
tional brush-type dc motors on an egual horsepower basis. 
Compared to its conventional counterpart, a brushless dc 
motor‘: structure is "inside out." In the brushless config- 
uration, the permanent magnets are attached to the rotcr and 
the ccenducting coils are placed in the stator. As its name 
Suggests, brushes have been eliminated from this tyre o£ 
NOCOT. Instead, ccmmutation of current in the stationary 
armature is accomplished electronically by switching cn the 
appropriate windings via transistors specially designed to 
provide the high currents needed for motor power. In crder 
to seguence the current to the coils, the commutation 


circuitry is provided rotor position feedback from sensing 


S 


devices, most commonly from Hall eftect sensere, These are 
simple semiconductor devices that develop a pclarized 
voltage depending upon the magnetic field passing through 
then. Small, highly sensitive and reliable, Hall effect 
senscrs reguire little power to operate and thus are widely 
used. 

Brushless dc motors possess all the advantages tradi- 
tionally associated with conventional dc motors such as 
linear torque-speed characteristics, excellent response 
times and high efficiency. In addition, there are a numcer 
of important advantages that brushless dc motors pcssess 
over standard brush-type dc motors. An obvious advantage is 
that electronic commutation means the elimination of comnu- 
tators, the wear and tear of brushes, and the build up of 
carbon dust which heretofore significantly limited the life 
of the standard dc motor. Since electronic current 
Switching essentially involves no component wear, there 
results improved reliability and ennanced lifespan. As weil 
as being subject to wear and tear, the action of brushes 
Sliding over commutators also generates radio freguency 
interference (RFI) which is a severe liabiiity in many 
applications. In fact, it is often in hostile and exrlosive 
environments that tke advantages of dc motors are most 
needed. Another benefit iS in heat transfer characteris- 
tics. By placing windings in the stator the thermal path to 
the environment is made shorter and more direct. with the 
removal of heat thus enhanced, mechanical and electrical 
malfunctions due to keat are reduced. A further important 
advantage of the brushless de motor is in the area of motor 
speed variability. Whereas ac motors usually operate at one 
speed fixed by the powerline frequency, any dc motor's speed 
can be ckanged to meet varying power requirements by sinply 
adjusting the dc input voltage. The electronic commutation 


scheme cf the brushless dc motor lends itself nicely to just 


10 


eiias type O£ variable speed control. eels: “a selats 


< 
(U 
}-4 

— 


eumgle Matter to attach to the Commutation circuit beard a 
small, inexpensive microprocessor whose logue can be 
programmed to run the motor at any of several speeds orc in 
either direction. A single motor can thus be proyranned, in 
place if desired, so that its speed, direction and torgue 
are matcted to a particular application. 

There are at least as many applications for brushiess 
dc motors as there are for brush-type motor systems as weil 
as a hest of new apfrlications. Of particular relevance to 


mimes thesis 1S the use of brushless dc motors in Ili 


Cy 
Er 
cr 


control actuator sytems of Navy tactical missiles, ania 
particular, the cruise missile. There are three general 
types of actuators that are in use in missile systems: pneu- 
matic, hydraulic and electromechanical. The maneuvering 
requirements of the mcre advanced tactical missiles call for 
high control torgues which heretofore could only be froducedi 
by high pressure pneumatic or hydraulic actuators. Bie awe 
these actuator systems are more susceptible to mechanical 
malfunction than their electrical counterparts. With the 
improved torgue to inertia ratios resulting from better rare 
earth magnet materials and the higher mecnanical reliability 
due to removal of brushes, as well as other advantages 
previously enumerated, the brushless dc motor is ina fosi- 
tion to compete with most pneumatic and hydraulic systems. 
This thesis is but one part of a detailed study of the state 
of the art in electromechanical actuators as applied tc the 
cruise missile system [Ref. 1]. 
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II. NATURE OF 


A. EROGRAMMNING LANGOAGE 


The objective of this thesis 


low-order mathematical 


the- response of a brushless de motor. 


Svstem Mcdeling Program (CSMP) was 


is a user-oriented ccmputer 


for modeling dynamic fhysical systems and is one 
It has the 


widely used in this area [Ref. 


response of 


Ze 
oc determining the 


€lther in a block diagram format or 


differential equations. 


Digital Simulation Language, 


model that will 


language specifically 


a systen 


While CSMP is an offshoot 


it embodies 


THE PROBLEM 


is to create a sinplifiied 


accurately simulate 


The IBM Continuous 


chosen for use since it 


created 
of the mest 
flexibility 
modeled 


that is 


as a set of 


ordinary 
of ISR: 


Fortran aS Des 


source language and retains much of Fortran's capability and 


CSMP is 
with 


£Vex aca eye 


intended 


an 
for use the 
utility and ease 

Simplified program ccntrol statements 
describe the mathematical equations or 


the system and from the flexibility of 


applications -orienmuca 
IBM mainfrane 


of using this program 


pregran 

The 
language stems fron 
tiated mos t 


systems. 


exactly 
physical variables of 


1ts program structure 


which contains preprogrammed function blocks that obviate 
the need for complicated Subroutine programming. In 
essence, CSMP was chosen because it permits the user to 
concentrate on the details of the physical system rather 


than on the time-consuming complexities of programming. 


Ee SYSTEM BLOCK DIAGEAM 


The block 
the CSMP progran 


diagram format was 


the input-output relationships cf a 


12 


chosen as tne 


Since it is a convenient 


basis for 


way of analyzing 


physical system and the 


flow of signals through it without having to refer to the 
system eguations. The basic idea of a block diagram stems 
from the application of Laplace transforms to a systen's 
integro-differential cr differential equations. IA eSet sect, 
a differential eguaticn is rendered into a simpler algetraic 
expression which, in turn, becomes a fundamental eguaticn o£ 
a system's functional blocks. Each block indicates ine 
‘relationship between its input driving signal and output 
response; ie., the transfer function of the input and 
uc out. The transfer function is defined as the ratio of 
the Laplace transforn of the output to the transform of tne 
input with all initial conditions assumed to be zerc. A 
Elock's transfer function does not include any information 
apout the internal structure of that part of a systen, 
merely the transformation of a signal between two pcints in 
a systen. Furthermcre, it should be borne in mind that all 
physical systems have certain transfer characteristics that 
are actually nonlinear to some extent. A de motor systen 
can be fairly accurately described using the transfer func- 
tion approach to modeling and, when nonlinear transfer char- 
acteristics exist, CSMP modeling can usually be obtained 


from preprogrammed ncnlinear function routines. 


C. DEVELOPMENT OF MCTIOR SYSTEM EQUATIONS 


In crder to construct the block diagram of a systen, 
one apfroach is to identify the various mathematical reia- 
tionships between the components of the system and to write 
balance equations that will identify each individual block 
comprising the system. The system equations for a brushless 
dc motor are presented in chapter three where all the 
deballs Of the Logie control and switching of transistors 
for ccmmutation of the motor are developed. What follows is 


a Simplified input-output characteristics model of the type 


le 


that has been a standard for brush-type dc motors used in 
control system studies [Ref. 3]. The development of this 
input-output model is the first learning step in the 
modeling process for the brushless de motor and has the 
additional usefulness of being a tool for a concurrent 
mechanical engineering study cf load toryue reyuirements for 
the MOtoE [Ree a): 

To write the electrical balance equation describing a 
kasze de noon, the classical loop method based on 
Kirchoff's voltage law is applied to the equivalent motor 


CcCiIfFEenit {Figurew2. 1) 





ed 
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Figure 2.1 Equivalent DC Motor Grrewit 
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The kasic dec notor has the coil windings on the armature and 
the magnetic field scurce on the stator. In the case of the 
brushless dc motor, the rotor consists of permanent nagnets 
that produce the magnetic field and the stator contains the 


Serlewindings. ithe Kirchoff Law for the stator circuit is 


es(t) = Ldi/dt + Ri(t) + eb(t) (eqn 2.1) 


where R and L are the stator resistance and inductance 
respectively, and eb(t) is the back enf. When a conductor 
moves ina magnetic field, or there is any relative mction 
ketween the magnetic field source and the conductor, a 
voltage is generated across the terminals of the conductor. 
iiche case of the de motor, the voltage is proportional to 
the shaft velocity and tends to oppose the current flow. 
The relationship between the back emf and the shaft velccity 


eS 


es (t) (Kn*) *wm(t) = Kb*wm (t) (eqn 2.2) 


es (t) Kb*dOm (t) /dt (eqn 23.3) 


EUBStltuting and rearranging equations 2.1 and 2.2 into 


their differential forms gives the following two equations. 


di/dt= (1/L) *es (t) - (R/L) *i (t)-(1/L) *Kb*wm (t) (eqn 2.4) 


don (t)vdt = wm(t) (eqn 2.5) 


The basic balance equation that describes the motor as 
ee Mechanical system derives £rom application of Newton's 
laws of motion to the system components. The dynamic eygua- 


tion for a motor coupled to its load is 


tm(t) = J*dwm(t)/fdt + B*¥wm(t) + tl(t) (egn 2.6) 


= 


lL 


where ti(t) is the load torque. ReabEbangings@tn]e cai47i one: 


becomes 
dwm (t) dt=(17J) *tm (€)— (175) tl (ef) Sea (eqn 2.7) 


Eis the total viscous friction coefficient of the mctor and 
is comprised of the sum of the viscous friction due to the 
motor, Bm, andthe load, Bl, as seen by the motor shaft. 
Likewise, J is the total system inertia comprised of the sua 
of the motor inertia, Jm, and of the load inertia, Jl, 
reflected through the coupling device to tne motor shaft. 
Foth Bl and Jl are functions of the motor's speed relucing 
mechanism where N is the speed reduction ratio (N>1). The 


following equations summarize the foregoing. 


B = Bm + Bl (eqn 227 
J = Jn + Jl (eqn 2a 
Bl = Bip/N {egn 2. 10) 
Jl = J1lp/N (eqn 2.11) 


Jlp and Plp are the inertia and viscous friction seen at the 
load side of the systen. Since a dc motor is essentialiyvea 
torgue transducer that converts electricai energy into 
mechanical energy, the following equation is necessary to 
show the relationshif amoung the developed torque, ta(t), 


the air gap flux, 0, and the Stator Gusment mite 


tm(t) = (Km*) *i(t) (een 2.12) 


Because the magnetic field in the motor is assumed uniforn 


and constant, this can be simplified to 


ile 


emtey = Kt*i (t) (egn 


KO 


. 13) 


where Kt is the torgue constant of the motor. 

PoUderwOnS waza te ChhOougiT 2.13 cepresent the cause and 
effect equations of the systen. MicmapplTecAatnon Sor eS (L} 
produces a current flow which causes the torcue and the back 
emf tc be generated. The torgue produced then causes the 
angular displacement Om(t). Given the differentiai equa- 
tions of the system, the Laplace transform can te aprlied 
and the block diagram of the system generated. As an 


example, the block fcr eguation 2.1 is generatec as follows 


es (t) - eb(t) = L*di/dt + R*i(t) (egn 2. 14) 
Aies(t) - eb(t)} =L(L*divdt + R*i(t)} (eqn 2.15) 
En(s) = Es(s) - EbD(s) = (LS + R)*TI(s) (Camnez =) 10) 


Mis)/7En(s) = Is/(Es - Eb) = 1t/(Ls + R) (eqn 2.17) 


Performing the same operation to the remaining differential 
equations, the block diagram in Fiyure 2.2 resuits. 

ities portant. to note that while a de motor is basi- 
cally an open-loop system, the back emf of this tyre of dc 
Motor acts as a natural feedback loop that tends to improve 
the stability of the sotor. 

Cnce the basic tlock diagram is developed, writing a 
Sole eeEogram aS simple and fairly straightforward. Care 
must be taken, however, to put the block transfer functions 
into the proper format for CSMP's functional blocks. A copy 
fetne  CSNP program for a basic dc motor is provided in 
Appendix A. 
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Prgure 2. 2 Block Diagram of DC Motor 
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Etiis CURVE EFREDEICT ION 


A. OVERVIEW 


Cne of the purreses of this model is to produce a set 
o£ performance curves that accurately simulate a given 
motor's behavior when operated under varying load cendi- 
mon S « As such the model can be a useful tool for studying 
the changing performance due to the configuration changes of 
a motor under develorment or for observing the behavior of a 
motor under varying conditions of application or environ- 
ment. Thus, having developed an input-output computer model 
for the DC motor, the next step was to assign values tc the 
model parameters and run the progran. Program inputs, 
constants and parameters were assigned values that are 
Mypacal of brushless DC motors commercialiy available. 
Table I provides the parameter values for a given motor that 
will ke used as a generic motor for the following analysis 


and discussion. 


TABLE I 


Typical Ccmmercial Motor Parameters 


Stator Resistance, R 2.74 ohms 
Seaton Inductance, L 0.0016 henries 
Horde Constant, kt 15.9 Oman ane 
pack ©§NF ~constant, Kb Qe 2 vyvolt/rad/s 
Rotor Inertia, Jn 0.001 OZ>1n/ se 
Per seOus FErerion 

Coerficient, En unknown oz-~in/srad/s 
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B. SPEED VS TGROUET CURVE 


Os 


Specifically, three curves Wereueprocice? S10 cue 
vS. generated torgue, current vs. torgue and output fower 
vs. torque. Since permanent magnet DC motors have linear 
characteristics, a straight-line curve resulted under ail 
load conditions when the parameters of Table I were inserted 
into the model and a voltage of 30 VDC was applied. In 
order to fully understand the individual contributions se. 
the varicus motor parameters on the speed-torgue curve, it 
was necessary to determine the mathematical relaticnships 
between the parameters. Referring to equations 2.1 and 2.5, 
the electrical and motor eguations when considered under 


steady state conditions become 


es (t) P*¥i(t) + Kb*wm(t) (eqn 24) 


tm (t) 


i 


Kt*i (t) (eqn 3.23 
Solving for i(t) and combining the equations resuits in 
es(t) = (R/Kt) *tm(t) + Kb*wm (t) (eqn 3.3) 


When there is no torgue applied to the motor, the current is 


small encugh to neglect and the no-load velocity becomes 
wnl = es/Kb (egn 3.4) 


Given that the applied voltage is a constant, it is apparent 
that knowledge of tke value of the back emf constant, Kb, 
Will previde the various no-load speeds for any value of 
applied voltage. On the other hand, in the absence of any 
known motor specifications, the no-load speed of a given 
motor can be empirically measured and then the back enf 


constant, Kb, can be mathematically derived. 
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Looking next at tke slope of the speed-torgue curve, it 
can be seen that whereas one endpoint of the curve termi- 
nates at no-load sveed, the other ends at stall which is 
that value of torgue which is large enough to literally stop 
the motor. From equations 3.1 and 3.2, with wm set to zero 


Bie, generated torgue, tm(t), at stall is ts given by 


ts = (Kt/R)*es (egn 3.5) 


Combining equations 3.4 and 3.5, the equation of a straight 


line for the speed-tecrque curve results in 


wo(t) = wnl - (Ro*tm(t)) {(egn 3.6) 
where 

Rm = R/(KD*¥Kt) (egn 3.7) 
and Km is the slope factor for the curve, also called the 
speed regulation constant. GbvVvousiy, Changing any or all 


of the parameters wiil alter the slope factor, Rn. HOWeVEL, 
Phe back emf constant, Kb, and the torgue constant, Kt, are 
strictly related by the common factor, aixc gap flux, $. In 
fact, when both constants are put inthe MKS system of 
units, Kb in volt/rad/s and Kt in Nm/amp, then they are 


equal as equation 3.8 shows. 


Kb(voltyrad/s) = Kt(Nm/amp) (egn 3.8) 


Thus, kncewledge of either constant provides knowledge of the 
other. Referring to equation 3.7, it can be seen that the 
final unknowns are the resistance, R, and the slope factor, 
Rm. If the resistance, R, is known from manufacturer speci- 


fication sheets, for example, then the slope, Rm, can be 
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easily sclved. On the other hand, if given an empirically 
measured speed-torgue curve, it 1S now obvious that the 
curve not only provides Kb which 18 Europ rovides Kty > laa 
also can provide a derived value for motor resistance, R. 
Cnce acquired, the three basic parameters, Kb, Kt and &, 
were plugged into the CSMP model and the speed-torgue curve 
for-a 30 VDC input was generated (see Figure 3.1). 

This same procedure can be used to produce a family of 
speed-torgue curves if it is desired to study a notor's 
performance under different input voltages. Figure (ere 
shows the family of curves produced when the values given in 


Table I were entered into the nodel. 


C. CURKENT VS TORQUE CURVE 


With the speed-tcrgque curve understood, the next task 
was to repeat the Same general procecure for the motor 
current vs. torque. As a starting point, the equation for 


the steady-state current 1s given Dy 
iss(t) = (es(t) - Kb*wm(t))/K (eqn 3.9) 


Since Kk and R were fixed by the speed-torgue curve, geéner- 
ating a current-torgue curve was straightforward procedure 


for any applied voltage (see Figure 3.3). 


D. CUTPUT POWER VS TCRQUE CURVE 


The final perfcermance curve to be studied was the 


output pcwer vs. torgue curve. The basic eguation fcr fower 


1s given by 


P(t) = tm(t)*wm(t) (egn 3.10) 


P(t) 


Kt*i(t) *wm (t) (eqn 3.11) 
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Figure 3.1 Motor Speed vs Torque Curve 


where pcwer iS measured in watts. Since this equation 
involves no new constants, fproducinjy the power curve w#as 
only a matter of adding a line to the existing nodeéel which 
wecorlngeuded a Conversion factor for changing the torgue 


units frem oz-in/amp to Nm/amp in order for the power to be 
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Figure, 322 Family of Speed-Torque Curves 
dimensionally correct in watts. Again Funning the pruejran 


with the parameters given in Table I, the model produced the 


power curve given in Figure 3.4 
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Figure 3.3 Motor Current vs Torgue Curve 


As expected both the speed-toryue and current-toryue 
Sumves Were linear. ithe power-torqgue Curve is obviously not 
Jinear, especially at high lcadis where the power begins to 
moeecifi fairly sharfly. Fron the modelling pcint of view 
this is due solely te the fact that the motor sicws decwn 


faster than the load torque increases (see ecuation 3.9). 
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Figure 3.4 Output Power vs Torque Curve 


Another factor that further accentuates this power rollcoff 
in an actual motor is the effect due to armature reaction. 
Whenever a current flows through the motor armature ina 
permanent magnet dc mector, the armature Lecones an electro- 
Dagnet which produces a flux that tends to oppose the ftiux 


produced by the permanent magnets. This has the effect of 
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partially demagnetizing the permanent hagnets. i aes 
demagnetization is a reversitle effect, meaning that as 
current returns to zero the permanent magnets return to 
their full strength. However, as the current increases and 
armature reaction sets in, it has the effect of decreasing 
the torque constant, Kt (as well as Kb). As equation 3.9 
magucates, this also contributes to the output power rolicf£ 
at high loads. Though current rermanent magnet fcc mctors 
made cf rare earth magnets have high coercivity that resists 
armature reaction, at high levels of rated current a small 
amount of armature reaction still occurs. For the purroses 
of this thesis, the model is considered to accurately siztu- 
dJate motor behavior fcr loads up to the peak power load tor 
the given applied vcltage. Prove tnaterOdde Voward, at 1s 
assumed that armature reaction is likely to occur and result 


in nonlinear behavior which 183 not included in this rodel. 


EF. MCTOR REVERSAL 


The final step in the modeling procedure was to reverse 
the motor's directior and ensure that mirror images or the 
three curves resulted. Teedo this; . Lt ewdsS Clrst necessary 
to replace the positive input voltage with a negative count- 
empart. Next, in order to maintain model consistercy where 
the lcad tcrque, tl(t), 1s treated as a torcue that offoses 
the moticn of the mctor, it was necessary to modify the 
model so that the load torque was a pcsitive value. This 
consisted of creating a CSMP procedure block that identified 
the input voltage as either positive or negative and then 
treated the load torgue accordingly. Once program tugs were 
removed, the mnodel accurately simulated motor reversal with 
proper speed, current and power values for the appropriate 


Mead conditions. 


on 


A. SWLITGHING ACTEON 


The next major step in this thesis waS to convert the 
basic model of a standard brush-type dc motor into a three- 
phase, four-pole brushless dc motor systen. Figure 4.1 
Shows this system with the three-phase stator windings 
configured in a standard 'star' connection with each winding 


oriented 120 electrical degrees from tne others. 


Figure 4.1 Controller Configuration 


The six transistors are connected to the ends of each stator 
leg and through logic-controlled switching action frovide 
three-phase full-wave motor coentrol. Figure 4.2 indicates 
the switching logic sequence for counterclockwise rotation. 
To produce clockwise rotation, the seguence is reversed. 
In either case, a pair of transistors will be switched at 
the start of each 30 degree (mechanical) interval which 


causes current to simultaneously build up in one leg, flow 
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Figure 4.2 Switching Logic of a 3-Phase Brushless DC MSotor 


steadily in a second, and decay to zero in a third. Given 
proper sequencing, the net developed toryue ideally reaches 
Geecteady State valte that causes motor rotation at a 
constant speed in the desired direction. 

For windings configured ina star arrangement, it can 
be seen that conduction iS ccntinuous in one leg while 
commutation occurs in the other two legs. For example, as 
Q1 and Q5 are energized between 30 and 60 deyrees (mechan- 
ical), current flows down leg A and into leg B. At the Same 
time, the current through leg Cc due to the energy storage 
behavior of an inductor decays to zero through diode, D3 
(see Figure 4.1). In the next sequence, 60 to 90 degrees, 
C6 is energized and ¢€5 switched off. Current has reached a 
steady state flow thrcugh leg A, but now it flows down leg C 
amomt ne GRErent ium leg 8 decays to zero through D2. In the 
hext sequence, 60°to 120 degrees, Q2 iS Switched on and Q1 
Switched off. Current has reached steady state flow through 
jeg C but now current builds up and flows from leg B while 


the current in leg A decays through ley C. This same action 


2g 


repeats in subseguent commutation states. LOOking at roang 
Cin Figure 4.1, it Can be Scenery stie- mec lGn- ices 
always continuous in one leq whether flowing into or fron 
node D, Kirchoff's Current Law reguires that the net current 
flowing in the other two legs must equal the fiow in the 
continuovs leg. This permits the following important 
simplification which was needed later in the stodeling 
process. Namely, the three-legged stator can be approxi- 
Mated by just two windings serially conrigured. Orne leg 
will always accurately reflect steady state flow conditions 
while the other will be treated as if i1t were in steady? 
state because the build up in one leg is balanced cry decay 
in the other. It is recognized that this model is a first 
approximation and that future modeles wili reguire addii- 
tional refinements in order to adequately represent the 
effect of switching transients on the transistor ani diode 


elements. 


Fe HALL EFFECT SENSOR FEEDBACK 


As waS mentioned in the introduction, the switching 
action of the commutation circuitry iS based upcn xTotoE 
position feedback from Hall effect sensors. in crder oe 
represent electronic switching and postional feedback, as 
well as the three phase winding configuration, it was neces- 
Sacry to modify the Flock diagram of Figure 2.1 to that in 
Figure 4.3 

Comparing Figures 2.1 and 4.3, it is seen that the 
chief difference is in the switching logic block and the 
addition of two more transfer function blocks to account for 
the additional windirgs. To better understand the logic 
relationship between the sensors and the current switches, 
the truth table for both counterclockwise and clockwise 


rotation iS presented in Table II. 
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EBegure 4.3 Equivalent Circuit With Electronic Commutation. 


A brief review of the action af a Hall effect senscr is 
dn order at this point. Consider a conducting material with 
a constant current-flow. When there is no ex tcr Nake mag ne tie 
field present, the current flows undisturbed in a straight 
line. But when a magnetic field is applied, a deflecting 


force, the Lorentz force given by Il x B (where I is the 
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TABLE a 


Sensor and Svwitcehingesoqac 


RES. RESwe RPS C Pa A PH 5B FE i 
1 1 Q high low open 
1 0 9 high open LOw 
1 0 1 open naserh icw 
Q 0 1 low high open 
0 1 1 low Open high 
0 1 0 open low high 

COUNTERCLOCKWIS& 
0 0 1 low high ofen 
1 0 1 open high l cw 
1 0) 0 high open Low 
1 1 0 high Low oren 
0 1 Q open low high 
0 1 1 Low open hier 
CLOGKWES® 


current, 1 is the length of the conductor, and 3 )2seu7 
magnetic field) causes the current to bend and electrons to 
pile up on cne side and positive charges, or holes, toa do 
the same on the other side. Thus a transverse Hall pcten- 
tial difference, Vh, develops. across the conductor (see 
Figure 4.4). 

Figure 4.5 shows that in the presence of a sinuscicacal 
Magnetic field produced by a rotating magnet, such as cne of 
the pole-pairs of a multiple-pole uotor, tke induced Hall 
voltage takes cn an ideal square wave shape. 

The logic circuitry of the sensor is configured in such 
away as to output a logic level 1 for positive Hall voit- 
ages and a 0 for negative voltages. Briefly reviewing the 
action of the motor: as the shaft turns, the rotor position 
sensors send this information to the switching logic -lock 
which then decides which phases are to be made high, lcw or 
open in accordance with Table II. Equivalently, the 


voltage, En, 1s then applied to the two appropriate windings 
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Figure 4.4 The Hall Effect 


and zero voltage 1s applied to the third winding in agree- 
Memt With the assumption that current is steady in one leg 
while the commutation in the other two legs is equivalent to 
a continuous current in one leg. As before, each voltage 


then produces a current. Because this 1S assumed to be a 
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linear system, the principle of superposition is applicable 
and, thus, the torgues developed by current flow through 
separate windings car be considered independently and then 
summed to give a total motor torgue. Beyond this point, the 
brushless motor compfenents are identical to its Erush-type 
counterpart and behave similarly. 
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ae OUDEL REVISION 


Ssipvewmodeling Of thisaversion of the brusniess dc motor 
called for two modifications to the brush-type nodel. 
Essentially, a CSME procedure function was constructed to 
act aS a subroutine that embodies the logic of the Switching 
aeti on. The logie cecognizes that aif the motor shaft is 
within scme degree rarge, for example 30 to 60 degrees, then 
it sets the appropriate rotor position sensors and energizes 
the switches for the appropriate two windings and deener- 
mgezes the switch for the third. In effect, only the forcing 
voltages are passed to the main program, whereas switch ane 


sensor positions are available for output (see Appendix A). 
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V. VARIABLE FLUX 


—_ Se ee ee = Se a 


A. AIR GAP FLUX 


The steady state behavior of both a brush-type and 2a 
krushless dc motor has been studied under the fundamental 
assumpticn that the krack emf generated is a product of tae 
shaft velocity times a constant, Kb (see Fguation 2.2). 
Likewise, the motor torque is equated to the product of the 
current times a constant, Kt (See Equation - 12) Bait ae 
Fquation 2.11 reveals in the case of the torque constant, 
Kt, the torque is in fact the product of some other constant 
cf properticnality, Kti, "times tie dina on. lu ore The 
Same is true for the tack emf constant, Kb; that is, kb is 
the product of some proportionality constant, Kb1, times the 
dll “qap flux, a0. In both cases, the flux, ¢, has beer 
assumed tc be constant which 1S in agreement with the 
general f~ractice when.modeling permanent magnet dc motors. 
In point of fact, the flux is not constant but rather is 
distributed sinusoidally, generally being the sum of a Ssiru- 
soid and one or more of its harmonics. How the motor is 
physically constructed determines this flux distribution. 
Two of the chief factors are the structure of the magnetic 
role faces and, most importantly, how the windings are 
distributed in the armature. The latter includes such 
factors as the number of slots, the number of coils fer slot 
and the spacing of the slots [Ref. 5]. The careful designer 
distritutes the windings so as to minimize the effects of 
harmonics. Regardless Of) tite ior constituents, these 
composite Sinusoids have an average value, of course, and so 
in practice it is this average value that becores the 
CONS tami line term, ¢. 


36 


Be. FLUX AS AN AVERAGE VALUE 


Pepoguce menuext s=stepeim the sdesign of a more realistic 
model was the determination of an average value for the 
Watible flux, avg, and the corresponding proctortionality 
Constant, KKi such that the ratio of dOavg and KX1 equals cne 
and the relationshifs Geeequattens 2.2 and 2a aig 
preserved. The eguations for back emf, eb(t), and develored 


torgue, tm(t), now become 


eb (t) Keo taewn.( Tt ) Gevanel 5): $9) 


moety == Ktp * 1{(t) (EGR Sc) 


where the back emf ccnstant, Ke pecdi@urio=cOrjue constant, 


Ktp, are now expressed as 
Kbp = Kb * (davg/Kk1) (een 528) 
Ktp = Kt * (@avg/KK1) (eqmn5 24) 


ASSUMing aSa first approximation that the air gap fiux 
varies aS a Simple sinusoid of unit magnitude, ee en a Se 
step toward finding the total system flux, Oavg, was to 
determine the relaticnship of the individual flux patterns 
in each of the three phase windings. Analysis of a four- 
fole magnet rotating under three windings separated by 120 


degrees (electrical) reveaied the following relaticnships 


da = sin(28 + 1/6) (eqn 5.5) 
db = sin(29 + 9176) (egn 5.6) 
dec = sin(26 + STY6) (eqn 5.7) 


ay 


where da, 65 and 0¢©€ are the Elux@§contrerou: cons or era 
respective phases. Feturning to an earlier assumption that 
the circuit can he approximated ads" tWom syindiggs 1:3 sera 
ignoring the third, then the average flux can alse be 
approximated from the sum of tne flux developed in two wind- 
ings at a time. This approximation iS further support cdma, 
the fact that the flux in the decay leg has, in fact, an 
average value of zero during every comuiutation interval. 
From the switching seyguence of Figure on, the two 
conducting legs were identified for each interval and, ‘£roa 
Eguations 5.5 through 5.7, their respective flux values 
calculated and summed. For every interval, the resulting 
fiux waveshape was identical and had an average value of 
davg = 1.631 (see Figure 5.1). 

Since KK1 merely normalizes davg, its value was easily 
determined to be KK1 = 0.613. The appropriate substitutions 


were then made in the CSMP model (See Appendix A for 


Revision Two). Since this amounted to nothing more than a 
rearrangement of constants, the model's output behavior 
remained unchanged. At this stage, the model simulates the 


behavior of a brushless de motor being electronically ccazau- 


tated and exhibiting constant, steady state developed torque 


th 


aie 


and motor speed due to a fixed, average value air gar 
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C. SINUSOIDAL FLUX 


The final stage of this thesis was to examine the 
Steady state behavior of the brushless de motor under vari- 
able air gap flux conditions. AS a fifSt approximaticn, sai 
flux was modeled as a Simple sinusoid and the flux relaticn- 
ships are those given in eguations 5.5 through 5.7 As 
before, only the flux relating to the two conducting legs 
waS considered. This composite variable flux, ¢v1, is shewn 
ID eigurer oc Equations 5.3 and 5.4 had to be modified to 


account for the different value of Elux 
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FLUX VARIATION 


2.0 


1.8 


1.4 


1.6 


ive 


COMPOSITE FLUX (WEBERS) 


1.0 


0 30 60 90 120 150 180 210 240 270 300 330 360 
ANGULAR ROTATION (DEGREES) 


Figure 5.1 Composite Plux Variation for Two Windings 
Kbpp = Kb * (Qv1/KK2) (egn 5.8) 
Ktpp = Kt * (0v1/KK2) (eqn 5.9) 


Recall tkat the objective was to keep the values of Kbpr and 
Ktpp as near as possible to those of Kb and Kt, respec 
tively. In order to roughly normalize Qvl, the constant 
value was determined as KK2 = 1.655. Since ?v1 was now a 
variable quantity, it should be recognized that Kbpr and 
Ktpp were no longer strictly constants and were better 


expressed as the follcwing functions. 


Kbpp = Kb (0v1) (eqn 5. 10) 
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Ktppr = Kt (Ov1) (eat oe td) 


The equations for back enf, eb(t), and developed torque, 


tm(t), now became 


eb (t, v1) Kbpp * wm(t) (eqn 5S. 12) 


tm(t,dv1) Kops melee) (eqn 5S. 123) 

Fecause the total flux was no longer a constant, 
average value, it was expected that variation in developed 
torgue (and, consequently, in motor speed) would result. 
Furthermcere, as the system bleck diagram of Figure 2.1 indi- 
cates, torque is ccnverted by the motor transfer function 
into a retational speed. This is a process that 1s equiva 
lent to inertia filtering of the torque ripple that resuits 
from variable flux. For the purposes of this thesis, riprle 
is defined as the ratio of the amount of variation from the 
maximum to the maximum itself. Having mace the apfrofriate 
Substitutions in the model (see Appendix A for kevision 
Three), the program was run under no-load conditions. Bie 
was crserved that thkere was, in fact, ripple in roth the 
develcped torgue and motor speed and that some inertia 
fii tering Wad oGeur. The next logical step was to icac the 
motor and study the effect on torgue and speed ripple. 
Table IiI indicates that under no-load conditions where 
speed was highest, the percentage of ripple in the output 
speed was at a minimum but that torgue ripple was at a 
maximum. 

It can be seen that as the load increased the amount of 
ripple in the develored torque decreased from about 84% to 
8% A closer exXamination of the torgue behavior uncer 
varying loads revealed that, even as the motor torgue necées- 


Sarily grew with increasing loads, in every case there was 
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PARLE = ie 
Torque and Speed Ripple Due to Sinusoidal Plux 


MOTCE . MOTOR fe ee Baa 

One COz = Lf) Ss Dh ey TORQUE Sf fe) 
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@eproximately a constant 8-9 oZ-in variation. This varia- 


meen Can be attributed primarily to the fiux variation in 
Ktp» and to a lesser degree to the variation in i(t) causec 
Bye the variation in Eack enf. Thus the magnitude of torgue 
variation remained akcut the same for all loads though the 
percentage ripple arrears to indicate otherwise. In the 
case cf motor speed, the amount of ripple apreared to 
increase slightly with increasing load and thus decreasing 
speed. In this instance, the magnitude of variation did 
increase a small amount with decreasinj speed - from 28 [pn 
fmemuo-to0aq to 52 rpm at 96 ozZ-ins - and suggests that speed 
ripple, at least, is somewhat dependent on load. Mis es 
due in part to the fact that as the motor sSiows, there is 
jess rotor momentum and thus less inertia filtering. In 
addition, at slower sreeds there are less frictional effects 
which at higher speeds also tend to filter out rippie. in 
general, it makes intuitive sense that at higher sreeds 
freyguency variations are harder to discern than at slower 


speeds. 


D. HARMCNIC FLOX 


Having observed the effects of Simple sinusoidal flux 
variation, the next step was to replace the single sinusoid 
With a flux composed of the sum of a sinusoid and its 


harmonics which is mcre representative of the air gap flux 
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Patterns in” an acral econ The back ems Of ~a EYL Gam 
commercial brushless dc motor is given in Figure 5.2 where 
the voltage between two terminalS was measured with the 
motor rotating at 1200 rpm 

A harmonics analysis of this waveshape shows that the 
principle harmonic was the fifth and that the waveshape was 


closely described by the expression 
3.0 Ssin{s): + “Og5o ssa nie) (egn 5. 14) 
When sukstituted into equations 5.25 thru 527, and the 


different angular speed taken into account, the follcyamg 


relationshirs for flux result 


da = 3.0*sin (20+ 1776) + 0.59*sin (100+5 1/6) (egn 5S. 15) 
db = 3.0¥*sin (29+9 1776) + 0.59¥*sin(106+9 7/6) (egn 5. 16) 
Oc = 3.0*sin(20+57T7E) + 0.59*sin (108+ 7/6) (egn 5.17) 


The same procedures used in the simple sinusoid model 
were then applied to the harmonics case wnere the compcsite 
flux changes to v2, a flux factor that varied to a slightly 
greater degree than did v1. Equations 5.8 and 5.9 haggage 


tke mcdified as follows 


Kbppp = Kb * (Qv2/XKK3) (egn 5.13) 


Ktppp = Kt * (OV2/KK3) {eqn 5.19) 
where KK3 approximately normalizes the flux and had the 
value KK3 = 5.38. As before, Kbppp ard Ktppp varied with 
the flux and so the new equations for back emf and develored 
torgue became 
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Figure 5.2 Back EMF Waveform Due to Harmonic Flux 
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Kb (Ov2) * wm(t) (eqn 5.20) 


eb (t, dv2) 


eb(t,dv2) Kbppp * wm(t) (egn 5.21) 


43 


tn (t, v2) Kt (Ov2) * i(t) (eqn 5.22) 


tm (t, 0v2) Kt ppp 1 i) (oqnu5.2s) 
The appropriate substitutions were oncS more fade in eee 
model (see Appendix A for Revision Four) and the modei was 


run under various loads. As anticipated, the variable flux, 


¢v2, produced ripple in the developed torgue and motor 
Speed. Table IV shews the resultant values for various 
loads. 

TABLE IV 


Torque and Speed Ripple Due to Harmonic Flux 


MOTCR MOTOR % RIPPLE 
LOAD (oz-in) SPEED (rpm) TORQUE SPEED 
0.0 2557 94.6 2.4 
aD 6 2087 38.5 2a 
64.0 1617 ONS 5.6 
96.0 1147 Vie 9.8 


Comparing Tables III and IV, it can be seen that the values 
for torque and speed ripple are of similar magnitudes. The 
somewhat higher values for torgue are due to the fact that 
in the dv2 case there was a 16-21 O0z-in variation as offosed 
to an average 8-9 cz-in variation in tthe $v1 case. The 
speed ripple behavior was much like the ¢v1 case with 
proportionate increases resulting from the increased torgue 
Lipple. Thus, for koth torgue and speed ripple due to v2, 
the model performed as expected; that is, slightly higher 
values of ripple occurred which were clearly attributable to 
the increased variability of the flux pattern of ccmpcsite 
Sinusoidal harmonics. 

An important observation made during runs of the v1 


model came as a consequence of the motor speed variability, 
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particularly when tke motor was unloaded. SinGe: Norload 
Speed was a function of Kb as given in eguation 3.4%, the 
vacriakle Kbpp had the effect of increasing the upper-end 
speed beyond the 2557 rpm baseline given in Figure 3.1 This 
increased magnitude ccupled with the swing of speed values 
created an undesirable effect in the unloaded conditior. In 
effect, at certain points in time the additional speed 
produced back emf values that exceeded the input voltace 
which in turn resulted in negative current and developed 
torgue values. Though these were transient negative values 
and had small effect on the steady state operation of the 
motor, still their presence was obviously unrealistic. The 
Peuuricn tO returning the motor to all positive values of 
current and torque was to increase the fricticn filtering 
effect by adjustment of the viscous friction coefficiert, 
Bm. It will be recalled from Chapter Three that a value for 
the Bm of the commercial motor being modelled was not avail- 
able and that for the basic model a value for Bm was derived 
from the curve fitting process. Thus, it seemed reasonable 
and permissable to repeat the procedure here. 

Since it makes sense that increasing the friction within 
a system will slow the system down, in this case increasing 
the value of Bm did just that. By increasing Bn ftrcn 
PeoegiS to 0.045 oz-in/yrad/s, the upper-end of the motor 
speed was reduced to a value that eliminated the undesirable 
negative values. This caused the no-load speed to te 
lowered slightly belcw the desired value of 2557 rpms. 
Recalling that no-lcad speed is fixed by the back enf 
constant, Kb, this parameter was adjusted from 0.1120 to 
0.1089 vclt/rad/s to tring the back up to speed. Since this 
adjustment was less than 3% and the manufacturer specifica- 
tions allowed a 10% +/- measurement error margin, this 
adjustment to Kb was permissable. OLeGOUrse, a propor- 


tionate adjustment to the torgue constant, Kt, was also made 
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since Kt and Kb must maintain 4 CoOnStart Srelaticusnizs = 
explained in Chapter Three. The effect on motor speed 
ripple was neglible. Reviewing the Ov2 model at this roint 
revealed that similar behavior in the no-load state occurred 
but cn a larger scale due to the greater variability of $v2. 
A sinilar set of adjustments was necessary to eliminate the 
negative current and torgue values in the unloaded cendi- 
tion. It required an adjustment of Bm from 0.00015 to 0.045 
Ooz-inyrad/s in order to ensure that all valueS were [fosi- 
tive. in addition, the back emf constant, Kb, was adjusted 
to 0.1192 volt/rad/s and a proportionate adjustment to Kt 
was Made. Again, all adjustments were within tolerances. 

Of course, it was earlier established in Equation 32.7 
that changes to Kb and Kt would change the slope of the 
speed-torgue curve. Thus, in both cases the slope was 
returned to its original value by adjusting the motor ré¢sis- 
tance, Rk. In the case of Sinusoidal variable flux, & was 
reduced trom 2.74 to Z.70 ohms. For harmoric variable flux, 
R was increased from z.74 to 3.0 ohms. Again, both ad juwae. 
ments were within the 10% +/- allowable deviation. It was 
also cbhserved that since both Kb and Kt no ionger remained 
constant but rather varied within small ranges, the slope 
waS not constant and varied accordingly in both the gOv1 and 
¢v2 models with the greatest change occurring when Kb and Kt 
Simultaneously reached the minimum values of their respec- 
tive ranges. Bearing in mind that these deviations in slope 
are fairly transient, their average values are most ifocr- 
tant and are the affrroximate values from which the final 


adjustments to R were based. 
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VI. SUMMARY 


Ree REMARKS AND CONCLUSIONS 


CSMP language is a convenient tool for modeling a 
rrushless DC motot. Positional sensor feedback and 
exmtching logic being Zunctions of time, the action cf e€lec- 
tronic commutation can be nicely simulated in CSMP. The 
model can be fairly easily modified in order to study cther 
fomor Configurations. For example, the three-phase star 
configuration used in this thesis could easily be changed to 
a grounded neutral, point Din figure 4.1, and the sane 
types of analysis performed. In future studies where the 
complexity of the mcdel is increased in order to include 
more design detail, the advantage of CSMP will beconte even 
more apparent. 

In the course of simulating the effects of different 
airc gap flux variaticns, it appeared necessary to make tinor 
parameter adjustments. in order to ofirset the feedbacx 
effects cf torgue and speed ripple. Within the scope of 
this thesis, adjusting the steady state performance was 
reasonable and justifiable. In a larger context, however, 
various control systems, particularly torgue generation 
Gomtrcliers, are specifically designed to reduce tcrgue 
ripple and its effects. So, in at least one sense, these 
adjustments were somewhat artificial. To the extent that 
parameter adjustment is necessarv and desirable, it must be 
understood that this is not a simple procedure since adjust- 
ment of one parameter aimost always requires adjustment of 
one cr more others. These adjustments in turn necessitate 
readjustzent of the criginal parameter. The process is thus 


an iterative one. This model permits this kind of parameter 
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adjustment, but it is a long and tedious@ process etia soe 


ke better performed Ly some sort of optimization algorithnn. 


PF. RECCMMENDATIONS FCR FUTURE STUDY 


The simulation cf ripple effects 1S particulary rele- 
vant to the control engineer who must design the control 
circuitry and power ccnditioner to meet specific performance 
reguirenents. This points to several areas for future 
study. One area is the control of torque generation and the 
several principle configurations used in brushless (ie 
motors. Two commonly used methods are the Sinusoidal torgue 
generation principle and the trapezoidal torque scheme. 
[Ref. 5]. In addition to torque generation control, a 
second area for investigation iS power control of brushless 
PC motors. In this thesis, output power was unregulated and 
waS studied simply as a function of a constant inyfut 
voltage. In practice, however, power is produced anid 
controlled by varying the supply voltage. Jarious schemes 
exist fcr power conticl; amoung these are iinear transistoc 
CORP EC. and pulse-width of pulse-ftrequency COn ties. 
[Reis 5]. 

A very important area that this thesis did not address 
was that of the switching transients that occur as a conse- 
guence of commutation. During each switching interval, each 
Winding is an inductcr that stores energy which causes a 
Slgnificant amount of current to flow at the instant o£ 
Switening. This can have serious effects, especially 1£ 
breakdown conditions exist amoung the commutation transis- 
tors and other contrel circuit elements. In this same area, 
Switching transients also occur within the Gontro! cCircupes: 
itself. Thus, a thorough investigation and understanding of 
the behavior of power transistors, diodes and other control 
elements used with electronic commutation needs to be done. 
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PEs InNG OF MODEL PROGRAMS 


A. BASIC PROTOTYPE MCDEL 


The program in this section iS a basic program that 
Slmulates the action of a standard brush-type DC moter. ae 
is the prototype from which the brushless DC moter model ani 


its several revisions are derived. 
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Be. REVISION ONE 


This is the first revision of the Basile op0ret y ce. 2 
model simulates the switching logic of a brushless [C motor 
which is based on feedback from Hall effect sensors. ED 
addition, the windings are treated independently and their 
contributions to developed torgue are superposed. Since the 
Superfesition results in twice as much current flow asa 
Single lumped coil, the total current is halved to retain 


the same overall motor behavior as in the prototype. 
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Com PREVESION Tae 


In the following program, the Noten” Seainr § Gauge) eee 
treated as the average value of the sum of the flux acting 
in two windings as explained in Chapter Five. FOr Cchven. 
lence of study, both windings are again lumped into a single 
winding as inthe frototype. This facilitates a closer 


focussing on the effects of variable flux on motor Lehavicr. 
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ENDJC6 
/* 


De. REVISION THREE 


The purpose of this version oi the program is treat the 
emecap flux aS varying in Simple Sinusoidal fashion. The 
total system flux is treated as the algebraic sum of the 


flux developed in the two active windings. 
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/* 


Mae REVISZCN FOUR 


This version of the model treats tne flux as varving 
according to the sum cf a sinusoid of fundamental frequency 
and its fifth harmonic as explained in Chapter Five. the 
total flux 1s again approximated as the algebraic sur of the 


flux developed in twce windings at a time. 
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APPENDIX B 


SAMPLE OUTPUT 


The following is a sample of the type of information 
tne mcdei can provide. Pieen is pAme CUbaL Gase, the rourth 
revision of the basic prototype, the harmonic flux case, was 
used. For completeness, the motor was run at loads froz 0.90 
up through 96.0 oz-in, which was the range of lcads with 
which all analysis was performed. In addition to such motor 
variakles as speed, current and power, a sampling of the 
feedback from Hall effect sensors as well as the switching 


action of the power transistors is included. 
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